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Abstract 

We derive the free energies of both the closed heterotic, and the unoriented, open and closed, 
type I string ensembles, consistent with the thermal (Euclidean T-duality) transformations on 
the String/M Duality Web. A crucial role is played by a temperature dependent Wilson line 
wrapping Euclidean time, responsible for the spontaneous breaking of supersymmetry at finite 
temperature while eliminating thermal tachyons, and determined uniquely by thermal duality. 
Conversely, we can show that the absence of a Yang-Mills gauge sector precludes the possibility 
of an equilibrium type II canonical ensemble prior to the introduction of background Dbranes 
or fluxes. As a consistency check, we verify that our results for the string free energy always 
reproduce the T 10 growth expected in the low energy field theory limits while displaying a 
dramatically slower T 2 growth at temperatures above the string scale. We present both the 
low and high temperature expansions for the one-loop heterotic and type I string free energies, 
results which follow from an explicit term-by-term evaluation of the modular integrals in the 
string mass level expansion. 
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1 Introduction 



T-duality invariance, the result of interchanging small radius with large radius, R — > a'/R, is a 
spontaneously broken symmetry in String/M Theory: in other words, a T-duality transformation on 
an embedding target space coordinate will, in general, map a given background of String/M theory 
to a different background of the same theory. To be specific, the circle-compactified E 8 xE s heterotic 
oriented closed string theory is mapped under a T-duality to the circle-compactified Spin(32)/Z 2 
heterotic oriented closed string theory [26, 28, 29, 37, 38, 6]. The circle-compactified type IIA 
oriented closed string theory is mapped to the circle-compactified type IIB oriented closed string 
theory [27, 30]. And, finally, the unoriented type IB open and closed string theory is mapped 
under a T-duality transformation to the, rather unusual, type I' unoriented open and closed string 
theory: the integer-moded open string momentum modes associated with the compact coordinate 
are mapped to the integer-moded closed string winding modes in the T-dual type I' theory [30, 45, 1]. 

The results in this paper are a beautiful illustration of the significance of the Euclidean T-duality 
transformations linking the six different weakly coupled string theory limits of the String/M Duality 
web at finite temperature. 2 It is clear that a Wick rotation on X° maps the noncompact 5*0(9, 1) 
Lorentz invariant background of a given supersymmetric string theory to the corresponding 50(10) 
invariant background, with an embedding time coordinate of Euclidean signature. The Wick-rotated 
5O(10) invariant background arises naturally in any formulation of equilibrium string statistical 
mechanics in the canonical ensemble, the statistical ensemble characterized by fixed temperature and 
fixed spatial volume (j3, V). The Polyakov path integral over connected world-surfaces is formulated 
in a target spacetime of fixed spacetime volume. Thus, the one-loop vacuum functional in the 
5O(10) invariant background computes precisely the sum over connected one- loop vacuum graphs 
in the target space R 9 in the finite temperature vacuum at temperature T [2]. The appearance of a 
tachyonic mode in the string thermal spectrum is an indication that the worldsheet conformal field 
theory is no longer at a fixed point of the 2d Renormalization Group (RG): the tachyon indicates 
a relevant flow of the 2d RG. The question of significance is then as follows: does the relevant flow 
terminate in a new infrared fixed point? If so, the new fixed point determines the true thermal string 
vacuum. An equilibrium statistical mechanics of strings requires that this fixed point belong to a 

2 The idea of pursuing the significance of Euclidean T-duality transformations on the different string theory limits 
of the String/M Duality Web came to me after reading Polchinski's discussion of thermal duality in the closed bosonic 
string theory given in his textbook [1], an approach I first considered in [19]. The observation that the one-loop 
string free energy has a first principles derivation from the Polyakov path integral, therefore expressible as an integral 
over the fundamental domain of the modular group, is apparent in Polchinski's analysis of the closed bosonic string 
ensemble in [2]; I simply applied that observation to the six supersymmetric string theories. For the heterotic string, 
I discovered that some of the necessary inputs for the one-loop vacuum amplitude already existed in the literature 
[6, 38, 37], except that their relationship to finite temperature string theory had not been noted. My presentation 
of the type I and type I' strings in [19] was guided by the well-established isomorphism to the heterotic strings [46], 
motivated also by the implications for the physics of the low energy finite temperature gauge theory limit [47, 48]. 
Finally, the impossibility of a type 11 canonical ensemble became apparent to me in the follow-up work [21], where 
I also discovered the significance of the temperature dependent Wilson line in both heterotic and type I ensembles, 
the connection to conjectures for strong- weak coupling duality without supersymmetry [39], and the T 2 growth of 
the string free energy at high temperatures. The realization that the conventional lore of the Hagedorn transition is 
flawed came to me in [21, 20], and was explored further in subsequent works [22]. Thus, the presentation that follows 
in this paper has been determined by pedagogy, rather than the chronology of my own understanding. 
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fixed line parameterized by inverse temperature (5 [12]: the precise analog under Wick rotation of 
the line of fixed points parameterized by the radius of a compact spatial coordinate in the 5*0(9, 1) 
vacuum. 

Target spacetime supersymmetry, and its spontaneous breaking in the thermal vacuum along 
the line of fixed points parameterized by /3, introduces new features into this discussion. We 
must require compatibility with the expected properties of the low energy field theory limit where 
the contribution from massive string modes has been suppressed, namely, those of a 10D finite 
temperature supersymmetric gauge theory. We must also require consistency with string theoretic 
symmetries of both worldsheet, and target space, origin in the Wick rotated SO(10) invariant 
background. In particular, Euclidean T-duality transformations must link the thermal vacua of the 
six different supersymmetric string theories in pairs: heterotic E s xE s and Spin(32)/Z 2 , type IIA 
and type IIB, and type IB and type I'. Remarkably, we will find as a direct consequence of the 
T-duality transformations, that the tachyonic thermal instabilities arising in all previous attempts 
to formulate an equilibrium supersymmetric string statistical mechanics in the canonical ensemble 
for the heterotic and type I strings are simply absent. For the type II superstrings, in the absence 
of Dbranes or fluxes, we will find that there is no canonical ensemble: the thermal tachyon free 
background also has target space supersymmetry, as a consequence of modular invariance. 

Having clarified that the canonical ensemble of heterotic and type I strings is well-defined at 
all temperatures, including the temperature regime far above the string mass scale, a'" 1 / 2 , we will 
establish several new results in this paper. First, we show that in either ensemble the growth of 
the free energy, F(j3), at high temperatures far beyond the string scale is only as fast as in a 2d 
quantum field theory. Thus, as was conjectured with only limited intuition as far back as 1988 
by Atick and Witten for closed string theories [41], there is a dramatic reduction in the growth of 
the free energy at high temperatures. More recently [1], Polchinski has shown that the T 2 growth 
in the free energy at high temperatures is a direct consequence of the thermal self-duality of the 
vacuum functional in the closed bosonic string theory. We will show in this paper that Polchinski's 
observation extends to both the heterotic closed, and type I unoriented open and closed string 
ensembles. In either case, the T 2 growth in the string free energy at high temperatures follows 
as a direct consequence of the Euclidean T-duality transformations that link the thermal ground 
states of the supersymmetric string theories in pairs. Furthermore, we will find that the full high 
temperature expansion in powers of f3 for the string free energy can be obtained explicitly upon 
term-by-term evaluation of the one-loop modular integrals in the string mass level expansion. 

Our starting point in this paper is the generating functional of connected one-loop vacuum string 
graphs, W(/3) = In Z(/3), where Z(/3) is the canonical partition function. W(/3) is derived from first 
principles in the Polyakov path integral formalism following [2]. The free energy, F(j3), and vacuum 
energy density, p(/3), can be directly inferred from W(f3). Let us recall the basic thermodynamic 
identities of the canonical ensemble [17]: 



Note that W(/3) is an intensive thermodynamic variable without explicit dependence on the spatial 
volume. F is the Helmholtz free energy of the ensemble of strings, U is the internal energy, and 
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p is the finite temperature effective potential, or vacuum energy density, at finite temperature. S 
and Cy are, respectively, the entropy and specific heat of the canonical ensemble. The pressure 
of the string ensemble simply equals the negative of the vacuum energy density, as is true for a 
cosmological constant, just as in an ideal fluid with negative pressure [17]. The enthalpy, H=U +PV, 
the Helmholtz free energy, F=U—TS, and the Gibbs function, also known as the Gibbs free energy, 
is G=U— TS+PV. As a result of these relations, all of the thermodynamic potentials of the string 
ensemble have been give a simple, first-principles, formulation in terms of the Polyakov path integral. 
Notice, in particular, that since P=—p, the one-loop contribution to the Gibbs free energy of the 
string ensemble vanishes identically! 

We should address the expected Jeans instability of a gravitating statistical ensemble [17, 15, 
41], and the basic definition of the thermodynamic limit in the presence of gravity. Consider an 
ensemble of total mass, M, and Schwarschild radius, Rs, with one-loop vacuum energy density, p 
= F^/R 0-1 . Recall that the Newtonian gravitational coupling, G N ~g^, where g s is the closed 
string coupling. Since M ~ pRg~ x , we have the relation: 



In other words, strictly speaking, it is only possible to take the infinite volume limit V — > oo, a' — > 
0, when the energy density per unit volume of the gravitating ensemble happens to be zero, or in 
the limit of weak string coupling [15, 41, 50]. 



2 The Impossibility of a Type II Canonical Ensemble 



We will begin by explaining why there can be no equilibrium type II superstring ensemble in the 
absence of Dbranes or background fluxes, either of which introduces Yang-Mills gauge fields in the 
low energy field theory limit of the finite temperature superstring theory. 

As explained at the outset in this paper, an equilibrium ensemble of type II strings requires that 
we identify a tachyon-free nonsupersymmetric background at all temperatures starting from zero. 
Let us examine some consequences of finding such a solution. Our interest is in a one-parameter 
family of nonsupersymmetric backgrounds of the type II superstring with Euclidean target spacetime 
R 9 xS 1 ; a line of fixed points of the worldsheet RG parametrized by the single parameter, T, the 
inverse of the circumference of the S 1 , whose physical interpretation is temperature. Continuity 
of the vacuum functional as a function of the parameter T requires that, at least for small values 
of T, W(T) take the form of an integral over the modular group of genus one Riemann surfaces 
with a modular invariant integrand. In addition, the leading terms in the closed string mass level 
expansion for small T must have a self-consistent interpretation in terms of the field theoretic modes 
of a 10D supergravity field theory at finite temperature. Thus, the target space supersymmetry of 
the zero temperature vacuum must be spontaneously broken at all temperatures different from 
zero. And the existence of a stable gravitating ensemble in thermodynamic equilibrium requires 
the absence of tachyonic thermal modes in the string mass spectrum. Finally, we must check that 
the vacuum energy density satistfies the criterion for the absence of the classical Jeans instability 
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in a gravitating ensemble [15, 41]. Thus, we can succinctly state the infrared consistency conditions 
required of an equilibrium ensemble of type II strings at finite temperatures: 

• the absence of thermal tachyons in the mass level expansion. 

• the spontaneous breaking of target space supersymmetry at low temperatures, T > 0. 

• the demonstration of a T 10 growth of the free energy at low temperatures, when we isolate 
the leading contribution to the string mass level expansion from the low energy, field theoretic 
super gravity modes alone. 

We will find that the first two conditions are incompatible: modular invariance of the one-loop type 
II vacuum amplitude turns out to be extremely restrictive. The spontaneous breaking of supersym- 
metry in the 10D type II superstring vacuum at temperatures different from zero necessarily implies 
the existence of a whole slew of low temperature tachyonic momentum modes. This invalidates the 
possibility of a type II string canonical ensemble, even at low temperatures far below the string 
scale. 

In order to understand the stringent limitations on the expression for the one-loop vacuum am- 
plitude in the finite temperature vacuum imposed by modular invariance, recall that the type II 
string mass level expansion results from a generic combination of the four, holomorphic and anti- 
holomorphic, Jacobi theta functions, weighted by, a priori, undetermined phases. In addition, the 
zero temperature mass formulae in each spin structure sector will now acquire additional contribu- 
tions from thermal winding and thermal momentum modes. Let us begin by examining the shifts 
in the ground state energy of the NS-NS vacuum of the type II superstring at finite temperature. 
Recall that this state is tachyonic in the zero temperature Hilbert space, eliminated from the mass 
level expansion for either supersymmetric type II superstring theory by a judicious choice of phases 
in the vacuum amplitude. It will be easy to verify that there is no corresponding choice of mod- 
ular invariant phases that can eliminate the tachyonic thermal modes of the NS-NS sector, while 
spontaneously breaking supersymmetry at all temperatures different from zero. 

Suppressing the oscillator contributions, the ground state energy in the NS-NS sector for physical 
states satisfying the level matching constraint takes the form: 



Notice that the pure momentum and pure winding states, (n, 0) and (0, w), are potential tachyons 
that enter into the general expression for the level expansion. In the absence of oscillator excitations, 
each pure momentum mode, (±n, 0), is tachyonic upto some critical temperature, T 2 = l/2n 2 n 2 a' , 
beyond which it turns marginal (massless). Conversely, each pure winding mode (0, ±u>), turns 
tachyonic beyond some critical temperature, T 2 = w 2 /8n 2 a'. Finally, recall that under a thermal 
duality transformation, the type IIA superstring is mapped to the type IIB superstring: 



thus interchanging the identification of momentum and winding modes in the type IIA and type 
IIB thermal specta. In other words, in the absence of a nontrivial Ramond-Ramond sector, the 
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result for the IIA and IIB one-loop vacuum functionals always coincides. In either case, upon 
compactification on the circle of radius 3/2ir, the zero mode spectrum takes the form [2, 1]: 
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so that the contribution to the path integral from the (n, w)th sector is: 
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Thus, upon compactifying either type II superstring on FfixS 1 , the one- loop vacuum amplitude 
takes the form: 
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Notice that the familiar choice of phase that projects out the tachyonic NS-NS vacuum at zero 
temperature, namely, the relative minus sign between (00) and (01) spin structure sectors, also 
suffices to eliminate all of the tachyonic momentum and winding modes from contributing to the 
string mass level expansion. This is the precise type II analog of the thermal tachyon spectrum 
of the closed bosonic string ensemble [2, 1], except that they have been rendered unphysical due 
to the choice of phase. In addition, we find that modular invariance has simultaneously forced the 
particular combination of worldsheet fermionic spin structures displayed in the expression above. 
Thus, target spacetime supersymmetry is unbroken! In other words, there is no viable type II string 
canonical ensemble. We need to introduce a new feature into the one-loop string vacuum amplitude 
that can loosen up the stringent constraints on the combination of fermionic spin structures arising 
from modular invariance. This new feature turns out to be the introduction of Yang-Mills gauge 
fields and, consequently, the possibility of spontaneous supersymmetry breaking via a temperature 
dependent Wilson line [19, 21, 22]. 



3 Free Energy of the Heterotic Canonical Ensemble 

Unlike the type II superstrings where the constraints from modular invariance proved far too re- 
strictive to permit a tachyon-free thermal ground state, the heterotic string theory has a Yang-Mills 
sector. This introduces the possibility of a Wilson line gauge background, loosening the con- 
straints from modular invariance while enabling a tachyon-free thermal vacuum and, consequently, 
a self-consistent formulation of the equilibrium canonical ensemble. Consider the ten-dimensional 
supersymmetric E 8 xE 8 theory at zero temperature. The a'— >0 low energy field theory limit is 
10D N=l supergravity coupled to E$xEs Yang-Mills gauge fields. What happens to the super- 
symmetric ground state of this theory at finite temperature? We wish to derive an analogue of 
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the zero temperature vacuum functional which describes a stable finite temperature ground state: 
thermal tachyons should be absent, target space supersymmetry must be spontaneously broken at 
finite temperature, and the free energy must grow as T 10 at low temperatures when only massless 
field theory modes are excited. Most importantly, since we wish to preserve the finiteness and 
perturbative renormalizability of the zero temperature ground state at finite temperature, it is im- 
portant to preserve the invariance of the one-loop vacuum functional under the modular group of 
the torus. Finally, we must require self-consistency with the thermal duality transformations: the 
E 8 xE 8 and Spin(32)/Z 2 heterotic string theories are related by the Euclidean timelike T-duality 
transformation: 

PesxEh -> $30(32) = 4ttV/73e 8 xE 8 , (™,w)e 8 xE 8 -> («' = W,w' = n) SO (32) , (8) 

thus interchanging the identification of momentum and winding modes in the E 8 xE 8 and SO (32) 
thermal spectra. Thus, the finite temperature vacuum functional is required to interpolate between 
the following two spacetime supersymmetric limits: in the /3— »oo limit we recover the vacuum 
functional of the supersymmetric E 8 xE 8 heterotic string, while in the (3=0 limit we must recover, 
instead, the vacuum functional of the supersymmetric Spin(32)/Z 2 heterotic string. 

It is helpful to begin by examining the resulting shifts in the ground state energy of the 
(NS+,NS+) vacuum due to the winding and momentum modes in the mass spectrum obtained 
upon compactification on the circle of radius fl/2ir. Suppressing the oscillator contributions, the 
ground state energy in this sector for physical states satisfying the level matching constraint takes 
the form: 
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where k^ belongs in the lattice E 8 xE 8 , and either E 8 lattice is spanned by the vectors: 
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(ni, • • • , n 8 ), (n± + |, • • • , n 8 + |), where all n; £ Z, and ^ ni £ 2Z . (10) 

i=l 

The corresponding result for the Spin(32)/Z 2 lattice takes the form: 

16 

(n 1: • • ■ , ra 16 ), (ni + |, • • • , rii e + |), where all n ; £ Z, and n ; £ 2Z . (11) 

i=l 

As for the type II superstrings, the pure momentum and pure winding states, (n, 0) and (0, w), 
are potential tachyons that enter into the general expression for the level expansion. In the ab- 
sence of oscillator excitations, each pure momentum mode, (±n, 0), is tachyonic upto some critical 
temperature, T 2 = l/2n 2 7r 2 o; / , beyond which it turns marginal (massless). Conversely, each pure 
winding mode (0,±w), turns tachyonic beyond some critical temperature, T 2 = w 2 /8ii 2 a'. Notice 
that retaining the relative minus sign between the contributions from (00) and (01) sectors of the 
right-moving superconformal field theory in the supersymmetric vacuum functional, as in the pre- 
vious section, would eliminate both 0((qq) l l 2 ) contributions to the level expansion, in addition to 
all of the (5 dependent tachyonic thermal modes. Thus, for example, all of the potential tachyons 
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in the (NS+,NS+) sector of the circle compactified supersymmetric ground state are eliminated 
from the level expansion for physical states in one fell swoop. But can this mechanism work in the 
absence of spacetime supersymmetry? We will need to identify a different chiral modular invariant, 
one that spontaneously breaks spacetime supersymmetry at finite temperature. The key lies in the 
introduction of a temperature dependent Wilson line background [19, 21, 22]. 

3.1 Axial Gauge and the Euclidean Timelike Wilson Line 



The generating functional of connected one-loop vacuum string graphs in the stable finite temper- 
ature vacuum will be given by an expression of the form: 



W het (P) =/5L 9 (4ttV)- 
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where the function Z\ iCt (f3) denotes the desired level expansion for the thermal mass spectrum 
of the heterotic ensemble at generic temperatures that can satisfy the interpolations to the zero 
temperature limit required by the thermal duality transformations. In particular, we wish to identify 
a suitable interpolating expression for Wh e t(P) which satisfies the infrared consistency conditions at 
generic values of (3, matching smoothly with the known vacuum functional of the supersymmetric 
E$xE 8 string theory at zero temperature ((3=oo): 
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while also recovering the vacuum functional of the T-dual supersymmetric Spin(32)/Z2 heterotic 
string in the presence of the T-dual Wilson line background: 
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It turns out that the desired expression for Zh e t(/3), and the pair of T-dual Wilson line backgrounds 
necessary for this interpolation, can be inferred using extant results in the heterotic string literature. 

The argument proceeds as follows. The modular invariant possibilities for the sum over spin 
structures in the lOd heterotic string have been classified, both by free fermion and by orbifold 
techniques [3, 8, 7, 9], and there is a unique nonsupersymmetric and tachyon-free ground state with 
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gauge symmetry SO(16)xSO(16). In the fermionic formulation, the GSO projection differs from 
that of the E s xE 8 heterotic string as follows: 



e irF+a+a' _ ■ 



-KiF+a'+a niF'+a+a 



+ 1 



(15) 



where a, a', and a, equal, +1, or 0, respectively, in the Ramond, or Neveu-Schwarz, sector for each 
of the three 8-fermion blocks: right-moving, left-moving gauge, and left-moving gauge' [7, 1]. The 
"all-positive" choice of worldsheet chiralities in every sector of the supersymmetric E$xEg theory 
are therefore altered to the more complicated, and nonsupersymmetric, projection: 



This choice of twisted boundary conditions introduces discrete torsion in the, alternative, orbifold 
description of this theory [8, 1], raising the ground state energy in the sectors that previously 
provided the massless chiral fermions of the supersymmetric heterotic string. Most illuminating, 
perhaps, is to understand this theory as the result of taking the noncompact limit of either circle- 
compactified E$xEg or Spin(32)/Z 2 heterotic string theories, upon introduction of a Wilson line 
background that spontaneously breaks both spacetime supersymmetry, as well as breaking the 
nonabelian gauge symmetry to SO(16)xSO(16). This interpretation is clarified by invoking the 
interpolating vacuum functional approach [32, 38, 37, 6]. 

Begin with the one-loop vacuum functional of the circle compactified nonsupersymmetric and 
tachyon free SO(16)xSO(16) heterotic string theory [7], where the S 1 is a circle of radius j3/2n: 



It is easy to verify by inspection that the integrand in this expression is invariant under the one- 
loop modular group. The expression also satisfies some of the infrared consistency conditions for 
a viable thermal ground state outlined in section 2: there are no thermal tachyons. However, the 
noncompact f3 — > limit of this expression does not lead to a spontaneous restoration of target 
spacetime supersymmetry since we recover, instead, the nonsupersymmetric 0(16) x 0(16) 10D 
modular invariant. We must incorporated a temperature dependent Wilson line background in 
order to interpolate between 10D supersymmetric and 9D nonsupersymmetric ground states. 

The spontaneous restoration of spacetime supersymmetry in the noncompact limit within the 
larger family of one-parameter continuously connected 9D 50(16) x SO (16) heterotic string vacua 
with nontrivial Wilson line gauge backgrounds was discovered in independent numerical investiga- 
tions by Itoyama and Taylor [38] and Ginsparg and Vafa [5, 6, 37]. The plots in the latter paper 
indicate clearly the distinction between taking the noncompact limit of the expression in Eq. (17): 
where the vacuum energy density goes thru a minimum at the self-dual point without touching zero, 



(NS+, NS+, NS+), (NS-, NS-, R+), (NS-, R+, NS-), (R+, NS-, NS-), 
(NS+, R-, R-), (R-, NS+, R-), (R-, R-, NS+), (R+, R+, R+) . 
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and taking the noncompact limit in the presence of a Wilson line background. It is in the latter case 
that the target spacetime supersymmetry is spontaneously restored at T = 0. Additional massless 
gauge bosons appear as we approach the T = limit, such that we recover the full SO '(32) gauge 
symmetry [37]. Recall that the gauge lattice vector in the equivalent worldsheet bosonic descrip- 
tion of Wh e t(P) will belong in a (1, 17)-component self-dual lattice of Lorentzian signature. Generic 
points in the lattice can be reached by Lorentz boosts that preserve the Lorentzian self-duality 
property, a pre-requisite for modular invariance, and each describes a heterotic string vacuum that 
can be reached by continuous interpolation from the theory above. The precise form of a generic 
vector in the self-dual lattice I^ 1,17 ) is given by [5, 6, 1]: 

(p%p° L ,k L ) - (p%P%,K) = (^-k-A -^A .A ,p L -k.A -^A .A ,k i+W ( £) A ) , 
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In particular, it is illuminating to consider the two lattice boosts that, respectively, interpolate 
between the nonsupersymmetric and tachyon-free 9D SO (16) x SO (16) vacuum functional and the 
pair of 10D E 8 xE 8 , and £0(32), supersymmetric string limits. The result for the required lattice- 
boosts is due to Ginsparg [6]. Starting with the supersymmetric 10D Spin(32)/Z 2 string, compactify 
on a circle of radius (3/2n with Wilson line gauge background: /3A° = 27r(l 8 ; s ). The closed form 
expression for the interpolating one-loop vacuum functional describing the exactly marginal flow 
to the 10D Spin(32)/Z 2 heterotic string vacuum can be found in the paper of Itoyama and Taylor 
[38] 3 : 



Z het (P) = ^(r)mr 8 [v(r)r W [ri(f)]- 4 x {(^o + S 1/2 ) 



+ (a, - 01/2) [|e* (e 2 6 + e^ 6 + ef) - ejelet - e^e 8 
+ {So - s 1/2 ) [-ie* (e 2 6 + e^ 6 + el 6 ) - e^e 8 e 8 + 

+ (0 O + 1/2 ) [-\Q\ (e 2 6 + e^ 6 + ef ) + e^e 8 e 8 - e^e 8 e 8 ] } . 



(20) 



As explained above, we use the T-dual of the standard basis for thermal mode summations in- 
troduced in [32], in terms of which the one-loop modular transformations are especially easy to 
verify: 



£ : n eZ, w even. £ 1/2 : n G Z + 1/2, w even. O : n G Z, w odd. 1/2 : n G Z + 1/2, w odd. (21) 



3 The vacuum functional of interest to us has been called the "Twist II" model by the authors of [38], as is clear 
from the plot in Figure 1 of their paper. The T-duality transform of the expression in Eq. (13) of Ref. [[38]] gives 
the desired expression for Zh c t(/3); a typo in the last line of Eq. (13) has been corrected. The flow to the 10D E$xE$ 
vacuum, instead, would be described by the interpolation with the T-dual Wilson line background, as in [6]. 
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A t^—1/t transformation leaves (£q+£i/2) and {Oq— C1/2) invariant, while interchanging {£q— £1/2) 
and {Oq-\-Oi/2). Under r— >r+l, the first three summations are invariant, while O1/2 maps to 
the negative of itself. It is easy to verify the invariance of this expression under a r — > — 1/r 
transformation. Under transformation, the first and third lines within curly brackets 

transform with a minus sign, while the second and fourth are interchanged with a minus sign. The 
overall minus sign is accounted for by the transformation of the eta functions. As before, w and 
n denote thermal windings, and thermal momenta, respectively. Notice the appearance of both 
half-integer and integer Matsubara frequencies along this flow. Notice, also, that because of the 
inclusion of a Wilson line background, the interpolating functional describing the flow to the T- 
dual E$xE 8 vacuum cannot be inferred directly from the expression above, because the T-duality 
acts nontrivially on the gauge field, in addition to interchanging momenta and windings. We must 
start afresh with the 10D E 8 xE 8 vacuum functional, compactifying on a circle of radius (3 with 
temperature dependent Wilson line background as before. The interpolation takes the form: 



+ (oo - O1/2) [|ej (e| + + e 8 ) 2 - e^e 8 - e'efe 
+ (£o - £1/2) [-§e* (el + e 8 + e 8 ) 2 - e^e 8 + eie 8 2 e 
+ (o„ + O1/2) \-W, (el + e 8 + e 8 ) 2 + e^e 8 - e*e«e 



} 



(22) 



It should be noted that in the presence of the timelike Wilson line backgrounds, the Yang-Mills 
gauge group is no longer a full SO (32), nor a full E 8 xE s , except in their respective noncompact 
limits: (3 — > 0. At generic values of (3, we have the gauge group SO(16)xSO(16)xU(l). Enhanced 
gauge symmetries of rank 17 will appear as we vary (f3,A (f3)), subject to the Lorentzian (17,1) 
self-duality conditions that ensure modular invariance. This is the usual phenomenon of Wilson line 
gauge symmetry breaking, a hallmark of heterotic string phenomenology. From the perspective of 
the low energy Yang-Mills gauge theory at finite temperature, the usual axial gauge quantization, 
A =0, has been replaced by a modified axial gauge quantization, which also spontaneously breaks 
the nonabelian gauge symmetry: / 9Ao=27r(l 8 ; 8 ) [19, 21]. We emphasize that the necessity for 
a Wilson line background in the finite temperature quantization arose as a consequence of our 
insistence on preserving the finiteness, and perturbative renormalizability, of the zero temperature 
ground state at finite temperatures. 

Let us derive from our expression for the string free energy the growth with temperature in the 
low energy field theory limit, namely, large (3. It will suffice to carry out this check in a single 
sector of the thermal spectrum, for instance, Co, with unit winding, but summing over all integer 
Matsubara modes, n G Z. We emphasize that modular invariance is not at issue here, since our 
interest is only in a field theoretic check of the scaling behavior as a function of T in the leading 
term of the string level expansion. Expanding the Jacobi theta functions in powers of qq = e~ 4nT2 , 
and extracting the coefficients of the 0(1) term in the level expansion of the zero thermal winding 
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number sector of Oq, gives the following leading contribution from massless target spacetime bosons: 



W)|(w=i ; n e z) = -2 • 2 8 • L 9 ( 47 rV)- 5 d n ■ (r 2 )" 




47r 2 a'n 2 (3 2 \ 
A-K 2 a' ) 



J 1 '* dnexp [2mnn] T (-5, \jl - r^j (23) 



The leading coefficient of the product of eta functions is {q l ^ 2 q) in heterotic string theory, and 
we isolate the compensating powes in the q expansion. Notice that all of the massless target 
space fermions of the zero temperature vacuum have acquired a tree-level mass that is linear in 
temperature, another self-consistency check for the low energy finite temperature gauge theory. 
The degeneracy at the massless level accounts for both the 64 bosonic states in the spin 2 gravity 
multiplet, as well as the 8.496 bosonic states in the spin 1 Yang-Mills multiplet. As an aside, as 
mentioned before in our discussion of the type II string level expansion in section 2, each degree of 
freedom in the spin 1 multiplet contributes with half as much weight as those in the spin 2 multiplet. 
The reason is that the level expansion counts target space degeneracies in terms of equivalent 2d free 
bosonic oscillator modes. The leading T 10 dependence in the free energy arises from the Matsubara 
spectrum; the subleading temperature dependence comes from the background of stringy winding 
modes. 



3.2 Thermal Duality and High Temperature Scaling Behavior 

The T 2 high temperature scaling behavior of the free energy can be inferred more elegantly by appli- 
cation of the T-duality transformation linking the thermal ground states of the two heterotic string 
theories. That argument proceeds as follows. In the closed bosonic string theory, the generating 
functional for connected one-loop vacuum string graphs is invariant under the thermal self-duality 
transformation: W(T) = W(T 2 /T), at the string scale, T c = l/2vra' 1 / 2 . It was pointed out by 
Polchinski [1] that we can infer the following thermal duality relation which holds for both the 
Helmholtz free energy, F(T) = —T ■ W(T), and the vacuum energy density, p(T) = —T ■ W(T)/V 
of the closed bosonic string ensemble: 

rp2 rp2 rp2 rp2 

F(T) = j^F(f), p(T) = ^p(f) . (24) 

In the case of the heterotic string, we will show that the thermal duality relation instead relates, 
respectively, the free energies of the E$xEg and Spin(32)/Z2 theories. Since we deal with a super- 
symmetric string theory, it is convenient to restrict ourselves to the contributions to the vacuum 
energy density from target space bosonic degrees of freedom alone: 

F(T) EsxEs = j^F( — )s P in(32)/Z 2 , P(T) Ef( xE 8 = ^p(^r)s P in(32)/Z 2 • (25) 
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Consider the high temperature limit of this expression: 

rp2 rp2 rp2 rp2 rp2 

lim p{T) E8xE8 = lim -j p(^)s P in(32)/Z 2 = Ijm ^p(^)spm(32)/Z 2 = 7^2P(0)spin(32)/Z 2 , 

(26) 

where p(0) is the contribution to the cosmo logical constant, or vacuum energy density, at zero 
temperature from target space bosonic degrees of freedom alone. Note that it is finite. Thus, at 
high temperatures, the contribution to the free energy of either heterotic ensemble from target 
space bosonic degrees of freedom alone grows only as fast as T 2 . In other words, the growth in the 
number of target spacetime bosonic degrees of freedom at high temperature in the heterotic string 
ensemble is only as fast as in a two-dimensional field theory. This is significantly slower than the 
T 10 growth of the high temperature degrees of freedom expected in the ten-dimensional low energy 
field theory. Notice that the thermal duality transformation interchanges the thermal winding and 
thermal momentum modes of the corresponding heterotic ensembles. This is the reason for the 
simple T 2 scaling behavior found in this duality relation. 

Thus, at high temperatures far above the string scale, the leading behavior is a significantly 
slower T 2 growth in the closed superstring free energy. Such a T 2 growth was first conjectured 
by Atick and Witten in 1989 [41]. Subsequently, it was shown to hold for the free energy of the 
closed bosonic string ensemble by Polchinski [1] as a direct consequence of its thermal self-duality 
property. 



3.3 High Temperature Expansion of Free Energy 



As a final bit of closed string methodology, we will now show how the modular integrals appearing 
in the result for the one-loop heterotic string vacuum amplitude can be carried out in closed form 
by the procedure of term-by-term integration. The result will be a power series expansion for the 
one-loop free energy of the canonical ensemble expressed as an exact function of temperature. 

We begin with a Poisson resummation on the infinite sum over winding modes, in order to put 
our expression for F(T)^ et in a form suitable for a high temperature expansion valid in the regime 
P — > 0. Denoting the coefficients in the level expansion by b^ et \ the result takes the form: 



-Fhet — 



dn f fW 4 E E E 

x &Pj ct) exp [-47rmr 2 ] exp 



/4vrVn 2 w 2 (3 2 ' 
-1/2' 



(3 2 



v ; 1 7-i/2 7 v /it^4t 2 2 2 UttV 



X 



6^ et) exp [-47rmr 2 ] exp 



'47rVn 2 ^ 
~J 2 



-7TT 2 



7r(p — UTi) 



+ 2ninwTi 



E E E 

=—00 p=— 00 m=0 

2 Aii 2 a r 



t 2 {3 2 



(27) 



Notice that the three integer- valued infinite summations in the expression for i 7 h e t(/9) given here 
account for, respectively, closed string mass level, thermal windings, and thermal momenta. It is 
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convenient to split the range of integration for the modular variable r 2 as follows: 



poo roc ry 

/ . dr 2 = dr 2 - 

J^j\-tI Jo Jo 



dr- 



2 , 



(28) 



where in the second term we must take care to express the integrand in a form that is manifestly 
bounded within the specified domain of integration. We will find that the former can be recognized 
as the Bessel function, K„(y), with a power series representation in the argument y: 



f 

Jo 



dx x y - 1 e--^ x 



v/2 



2 1-1 K v ( 2^7 ) , Re 5 > 0, Re 7 > 



where 6 ~ ' ' < - — - - ' 1' 2 



P 2 



-(p-nn) 2 , 7 = 47rm + 47rVn7/r, and u = -11/2 . (29) 



The latter can be expressed in terms of the standard Whittaker functions, W_*+i k(u), with integral 
representation: 



[ U dx x"- 1 e- 8 ' x = S^-^u^^e-^W^ „ (5/u) 

JO 2 >2 



47r 3 o/ - 1 1 i 

where 5 = (p — nrx) 2 , u= y 1 — r 2 , and v = h k, k G Z + , (30) 
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following substitution of the Taylor expansion of the function exp[7T2] in the integrand. Note that, 
with this substitution, we can verify that the integrand is always bounded in the region \r 2 \ < 1, 
term-by-term in the Taylor expansion. Substituting in the expression for the free energy given in 
Eq. (27) we obtain the result: 



-Fhet — 



fill POO r\ Tn OO OO OO 

L»(4.V)-/ in £ £ E 

J -1/2 •>V 1_T l 1T 2 n=-oo w=-oo m=0 



x exp [-47rmr 2 ] exp 

,1/2 



4n 2 a'n 2 w 2 (3 2 



-7TT 2 



(3 2 



/L/z roo 
d n 
-1/2 Jy/l-T? 



x &J^f exp [-47rmr 2 ] exp 



! rl/2 00 00 00 

^tV)" 6 • /T 1 / dr x £ £ £ 6 £et>(p _ nTl) -n/2 



-7TT 2 



* f wJ +2 ™ Tl . 

1/2 00 00 00 

E E E 

n=— 00 p=— 00 m=0 

An 2 a'n 2 \ . 2 4ir 2 a'~ 
* ir(p — nri) 



_ 4 / r 2 
?4r 9 2T2 V47r 2 a 



(3 2 J ^ LJ r 2 (3 2 
00 00 00 



n=— 00 p=— 00 m=0 



X 



2 , 

7rcr 



11/4 



8ttV 



11/2 



2 {p-nn) n H 



2\ V2 X 



I rl/2 00 00 00 00 

\L\4n 2 aT 5 • /T 1 / ^ E E 



n=— 00 p=— 00 m=0 fc=0 



L(het) 

~kT 



X 



4ttV 



n 2 + 



m(3 2 ^ 
ir 2 a' 



4ttV 



(3 2 



(p - nnf 



(fe-13/2)/2 



(1 _ r 2)(fc-9/2)/4 



13 



x exp 



2ttV (p -nnf 



(fc-9/2)/2,(fc-ll/2)/2 



An 3 a' (p - riTif 



(31) 



Expressing, in turn, the Bessel, Whittaker, and exponential functions as convergent power series 
expansions in their respective arguments, enables all of the T\ dependence in the integrand to be 
extracted in the form of an algebraic power series. The term-by-term integration over T\ can then be 
carried out explicitly. Note the useful relations W\^(y)=W\_^(y) , K u (y)=K_ v (y). We will make 
use of the following power series formulae: 



„ , v [T _ * (5 + r)! . , r 8ttV. . ( 2 m(3 2 \ 1/2 

K 11/2 (y) = ^Yy 6 r ? r!(5-r)! (%) ' and V = ~^ ~ nTl) ' n + 



P 2 



ir 2 a' 



.(32) 



Notice that, term-by-term, the T\ dependence in this series is remarkably simple. Also, notice that 
the second subscript, fi, in the Whittaker function is such that 2/i is half-integer. We can use a 
functional relation, and the power series representation: 



M 



r(-2//) 



»--&)■ 



M 



z + 



r(2//) 



r(| + /i-A) 



where : = (P ^ \ = -{k- 9/2)/2, /i = (k - ll/2)/2 



/? 2 



(33) 



Prior to r x integration, it is convenient to include the Taylor expansion of the exponential, e z , as 
well as the powers of binomials, (p — nri), (1 — i -2 ) 1 / 2 , accompanying the Whittaker function in the 
expression in Eq. (31), defining the coefficients, w^l^((3,a'), and numerical subscripts, o±(A,/i), as 
follows: 



(p - nn)( fc - 13 >/ 2 (l - r 2 )( fe - 9 / 2 )/ 4 e-2 z W_ {k _ 

r+s+a± T(— 2/i) ( + ) 

r(i-^_A)>' s ) 



(fc-9/2)/2,(fc-ll/2)/2l^J 



EE-r 

r=0 s=0 



r(2/i) 



r(i + / i-A) 



(-) 

(r,s) 



(34) 



A similar expansion holds for the term proportional to Kn/2(y), with the coefficient n(a', 13) defined 
by the relation, y = (p — nri)/c. The result for Fh e t consequently takes the form: 



1 oo oo oo 

ft* = -^L B (47rV)- B E E E 4 het) 

n=— oo 

/■1/2 

x / drAp — nr — 1) 

7-1/2 



n=— oo p=— oo m=0 
,11/2 



2 . 

n + 



2T11/4 



7r 2 a' 



7T 



2n(p — riTi) 



oo oo oo 



E E EE 



L(het) 
m 

~kT 



X 



1/2 
■1/2 



rfr n E E r i 



,r=0 s=0 



n=— oo p=— oo m=0 fc=0 

r (~ 2 ^) x ( + ) 



r=0 
2 , ™/3 

n + 



7r 2 a' 



(-) 

(r,s) 



r(2/x) 



r(i + / ,-A) 



(3fe-13/2)/2 



(35) 
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With these substitutions, all of the T\ integrals appearing in our expression for Fh e t can be evaluated 
explicitly. Recall that m and n denote heterotic string mass level, and thermal momentum number, 
respectively. The infinite summation over p resulted from the Poisson resummation of thermal 
winding numbers, necessary to cast the expression for the free energy in a form suitable for a 
high temperature expansion in (3. Thus, following the one-loop modular integrations, we obtain an 
exact high temperature expansion for the one-loop string free energy as a function of f3, and the 
degeneracies, b^ et \ in the heterotic string mass level expansion. 

Finally, it should be noted that, using the power series representation in j3 given above, one 
could demonstrate the analyticity of infinitely many thermodynamic potentials for the heterotic 
ensemble. Namely, potentials defined by taking arbitrarily large number of temperature derivatives 
in the vicinity of the self-dual temperature. This is clear evidence that the thermal duality transition 
linking the high temperature behavior of one heterotic ensemble to the low temperature behavior 
of the other ensemble lies within the Kosterlitz-Thouless universality class [18, 35, 34, 28, 20]. 

4 Free Energy of the Unoriented Type IB Ensemble 

The free energy of the unoriented open and closed type IB string ensemble at one-loop order in the 
string coupling receives contributions from surfaces of four different worldsheet topologies [1] : torus, 
annulus, Mobius strip, and Klein bottle. Given our analysis of the heterotic and type II superstring 
ensembles in previous sections, we expect that the Yang-Mills fields will play an essential role in 
enabling discussion of an equilibrium canonical ensemble. As explained in section 2, pure type 
II closed string theories have no Yang-Mills gauge fields, but upon coupling to Dbranes, the full 
unoriented open and closed string theory indeed contains timelike Wilson line backgrounds. This 
is the case we shall consider in this section. 

The torus contribution sums over closed oriented worldsheets, and the result therefore follows 
from the analysis of the type IIB superstring described in section 2. In particular, there are no 
tachyon-free thermal ground states that break target space supersymmetry and, consequently, no 
equilibrium canonical ensemble. The functional, W t0 r(P), therefore vanishes identically. Thus, at 
one-loop order in string perturbation theory, the oriented closed string sector will not contribute to 
the free energy of the type IB thermal vacuum. 

Thus, our focus must shift to the open and unoriented string sectors of the type IB thermal 
spectrum. An equilibrium description of type IB string statistical mechanics in the canonical en- 
semble requires a tachyon-free thermal spectrum in the full temperature range, in addition to the 
absence of massless Ramond-Ramond sector tadpoles. The vacuum functional at one-loop order in 
the presence of iV = 2 5 D9branes can be written in the general form: 



where the function Zjb(/3) denotes the level expansion of the type IB unoriented thermal ensemble, 
with contributions from annulus, Mobius strip, and Klein bottle. 




(36) 
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Notice that there is no constraint analogous to modular invariance on the open string mass 
spectrum. In addition, there are no thermal winding modes. Consider the contribution to the path 
integral from the nth Matsubara mode: 

e -4*Wt/p^ n, n + | G Z . (37) 

Focussing on the t— >oo IR limit of the string amplitude at first, at low temperatures (3 — > oo, we 
see that the field theoretic modes with low values of n dominate. But in the opposite t — > limit, 
we must include arbitrarily high values of n in the Matsubara spectrum. What about the high 
temperature limit, with (3 — > 0? Even in the t — > oo limit, modes with arbitrarily high n are now 
important. Thus, a thermal duality transformation to the Euclidean T-dual type I' ensemble would 
illuminate the high temperature behavior of the ensemble. As in the case of the heterotic string 
ensemble, the key will lie in introducing a temperature dependent Wilson line gauge background 
in order to achieve the spontaneous breaking of supersymmetry without the introduction of ther- 
mal tachyons. We must simultaneously require that Z(ib)(/3) preserve the Ramond-Ramond sector 
massless tadpole cancellations necessary for infrared finiteness and perturbative renormalizability 
of the finite temperature vacuum. 

Recall that the Wilson line introduced in the Spin(32)/Z 2 heterotic string breaks the nonabelian 
gauge group to 50(16) xSO(16), while spontaneously breaking supersymmetry. Let us compactify 
the type IB SO (32) string on a circle of radius (3/2ir, turning on the temperature dependent Wilson 
line [30, 45, 1]: 

(3A = (1 8 ,0 8 ), 9 l = 2n, i = 1, • • • , 8; ft = 0, i = 9, • • • , 16 . (38) 

This is the SO(16)xSO(16) type IB vacuum, and the Wilson line wrapping the X° coordinate has 
spontaneously broken supersymmetry in addition to a partial breaking of the nonabelian gauge 
symmetry in the supersymmetric SO (32) vacuum. Thus, we wish to find the finite temperature 
vacuum functional that interpolates between the tachyon-free and nonsupersymmetric type IB 
SO(16)xSO(16) string at finite (3, and the supersymmetric type IB 5*0(32) string in the limit 
where the Wilson loop is infinite radius. 

The Euclidean T-duality transformation, (3 — > (3'=4ir 2 a'/ (3, maps the configuration of D9branes 
coincident with an 09 plane to a configuration of a pair of 08 planes, each coincident with 16 
D8branes, at the two endpoints of an interval of length (3'. The thermal modes of the type IB 
ensemble correspond to a Matsubara-like frequency spectrum with timelike momentum: p n =2nn/ (3, 
where n, n + \ e Z; for the T-dual type I' ensemble this role is played by the infinite summations 
over thermal winding modes. Note that the D8branes are each paired with their images, so that 
the gauged background has a total of 8 2 oriented stretched strings of length f3' . Supersymmetry is 
spontaneously broken due to the presence of the stretched strings, which contribute with positive 
thermal tension to the vacuum energy of the finite temperature vacuum: 

T((3) = (3 2 /An 2 a' , (39) 

contributing an overall shift in the masses of the thermal tachyons in the NS sector. The GSO 
projections on worldsheet fermions eliminates all tachyons from the physical state spectrum, as in 
the supersymmetric vacuum. 
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The cancellation of both RR-RR, and NS-NS, sector massless tadpoles, and the absence of open 
string tachyons, in a nonsupersymmetric and tachyon-free 9D SO(16)xSO(16) type IIB orientifold, 
was shown by Blum and Dienes in [39]. As a consequence, both the cylinder and the Klein bottle 
amplitude vanish identically in the 9D vacuum, and the result for the vacuum energy density can 
be written entirely in terms of the contribution from the Mobius strip: 



yNS-NS 
^Mob 



7R-R 
J Mob 



(/*) = 




6qi (it; 0)010(^0)' 
r)(it)&oo(it) 

8oi (it; 0)010(^0)' 
rj(it)Qoo(it) 



(40) 



The authors of Ref. [39] have computed the individual contributions from the Mobius strip and 
torus amplitudes numerically down to values of the radius at which the divergence in the torus 
amplitude sets in, as a consequence of the closed string (high temperature) thermal tachyon. The 
asymptotic approach to the supersymmetric SO (32) vacuum at low temperatures is evident in their 
plot. As explained earlier, we believe it is more meaningful to interpret these results as evidence 
that the worldsheet RG flow is towards the IR stable supersymmetric vacuum in the closed string 
sector: the one-loop vacuum amplitude computes the tree level mass spectrum, and at this order in 
string perturbation theory the supergravity and Yang-Mills sectors have not as yet "communicated" . 
Thus, supersymmetry remains unbroken in the gravitational sector, and the nonvanishing one-loop 
vacuum energy density of the type I model comes wholly from the Mobius strip. 

We believe the correct physical interpretation of the Blum-Dienes result is as follows. Unlike 
the heterotic case, in the type I strings, there are thermal interpolations between the 10D super- 
symmetric SO (32) heterotic vacuum and either the 9D nonsupersymmetric and nontachyonic type 
IB, or the T-dual type I', SO(16)xSO(16) string vacua. The former has only thermal momentum 
modes, and the latter only thermal winding modes. The result for the free energy of the type I' 
ensemble takes the form: 



F V (P) = -L 9 (4ttV)- 



dt 
2t 



x 



E 



e i(it;0)eio(it;0) 
r](it)Q 00 (tt) 

a -2irf3 2 w 2 /iir 2 a' + e -2TT(3 2 {w+\) 2 }/^ 2 a' 



(41) 



or the Euclidean T-dual transformed result for the type IB ensemble: 



6 i(it;0)eio(it;0) 



1 4 



00 

x E 

n=— 00 



r}(it)Q 00 (it) 
+ e 



8n 3 a'n 2 /(3 2 j_ o -87rV(n+i) 2 ]//3 2 



(42) 



Thus, we have found evidence for a stable equilibrium ensemble of type I unoriented open and 
closed strings in the limit of weak string coupling. 



17 



4.1 The Low Energy Limit of F m (T) 



At low temperatures far below the string mass scale, P>>a n ^ 2 , we expect not to excite any thermal 
modes beyond the lowest-lying field theoretic modes in the open string spectrum. A useful check of 
self-consistency with the low energy finite temperature field theory limit is to verify the expected 
T 10 growth of the free energy. For simplicity, we will restrict ourselves to the sector of the thermal 
spectrum with half-integer moding of Matsubara frequencies. An analogous expression holds for 
the sector with integer moding. 

We begin with the one-loop contribution to the type IB string vacuum functional given in Eq. 
(42), including both integer and half-integer Matsubara frequencies. Expand the integrand of the 
modular integral in powers of q=e~ 2nt , thus isolating the t^oo asymptotics of the modular integral 
which is dominated by the massless modes in the open string spectrum. We have: 



e 01 (it;0)eio(rt;0) 
r)(it)Qoo(it) 



x 



E 



e -8^a'n 2 //3 2 + e _ 8 ^ a '(n+|) 2 //3 2 



(43) 



Notice that at sufficiently low temperatures, an infinity of Matsubara modes can contribute to the 
(3 — > oo limit of this expression since the integral on t ranges over [0, oo]. Thus, we must include 
the full sum over the full Matsubara spectrum when computing the leading low energy field theory 
limit of the string free energy. 

We recognize in our expression the integral representation of an Euler gamma function with 
negative argument, T(— 9/2), which may be defined by either the product formula for gamma func- 
tions, or by analytic continuation in the argument of the gamma function. Thus, performing the 
explicit integration over t for the leading term in the level expansion on to, and restricting to the 
summation over just the half-integer moded Matsubara spectrum for illustration, gives the result: 
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(44) 



where L 9 denotes the nine-dimensional spatial volume. T(—n) would have to be defined by analytic 
continuation, invoking a complex integral representation. The factors within the second pair of 
square brackets in Eq. (44) result from the integral over the world-sheet modulus, t. The modular 
integration is followed by the infinite summations on even and odd integers, giving the result within 
the first pair of square brackets. The summations have been expressed in terms of the Riemann 
zeta function ((z, q): 
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(45) 
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and B n {x) is a Bernoulli polynomial, B' n {x)=nB n -i{x). The parameter, (3q, has the dimensions 
of an inverse temperature, and it characterizes the asymptotic limit of the free energy for the low 
energy supersymmetric gauge theory. Recall that the string vacuum functional is dimensionless. 
Thus, we have demonstrated that the one-loop free energy, F(/3)——Wib(/3)//3, grows as T 10 in the 
asymptotic low temperature regime. This is precisely the behavior expected of a ten-dimensional, 
finite temperature supersymmetric gauge theory. Explicitly, we have: 

/?o = [2 4 ( J B 11 (|)/ll)r(-5)] 1 ^(47rV) 1 /2 . (46) 



4.2 High Temperature Expansion 



Conversely, when we include the contribution from all of the massive string modes, we find a T 2 
growth in the string free energy at temperatures far above the string scale. Set t—*l/t in the 
expression for the level expansion, and expand in powers of e^l 1 . The ultraviolet asymptotics 
of the mass level expansion dominates the high temperature regime, giving the following leading 
contribution to F m ((3): 
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(47) 



Thus, we find that the free energy at one-loop order in every perturbative string ensemble, whether 
open or closed, bosonic or supersymmetric, grows as T 2 at temperatures far above the string scale, 
a conjecture originally made by Atick and Witten [41] for closed string theories alone, and later 
shown to be a consequence of thermal self-duality in the closed bosonic string theory by Polchinski 

[!]■ 

Including all of the corrections to the ultraviolet asymptotic limit of the open string mass 
spectrum, we find that the term-by-term integration over the modulus, t, can be expressed in terms 
of the Bessel function K\(x). The free energy takes the form: 



n=0 J ° 2t 



dt 
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(48) 



where the are the coefficients in the mass level expansion for the type IB thermal mass spectrum, 
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obtained by expanding in powers of q=e 2w ^: 

Performing the integration over t results in the expression: 
F m {(3) = -L 9 (4n 2 aT\Pi/P) 2 

oo oo 



-2irm/t 



(49) 



m=l 



■ x K x {x) 



(50) 



where P 2 x 2 =64n 4: a' m(n + ^) 2 identifies the argument of the Bessel function. We emphasize that this 
result for the type IB free energy is valid at all temperatures, and includes all of the massive string 
modes, as well as the full half-integer moded Matsubara frequency spectrum. A similar expression 
can be derived for the integer moded frequency spectrum. The factor in square brackets in the first 
equality arises from the modular integration. The infinite summation over thermal mode numbers 
therefore yields ((— 1, |). The large x, or high temperature, asymptotics of the Bessel function 
indicates clearly that all of the corrections to the leading T 2 growth of the string free energy are 
exponentially suppressed. The asymptotic expression for the one-loop free energy therefore takes 
the form: 

F IB (T)/L 9 (47rV)- 5 = -[3 2 T 2 - £ £ A m T e -^ m ^ T . (51) 



m=l n=l 



where the leading T 2 growth is found to be corrected by a negligible Te T dependence, at high 
temperatures. Recall that Pi takes the form: 



(2 7 r)2 4 r(-l)fi 3 (|)/3l^ 



(52) 



where we have substituted for ((— 2, |). As before, r(— n) has to be defined by analytic continuation. 
Notice the minus sign in the coefficient, reversing the sign of the free energy, and clarifying that the 
dominant UV contribution is from target space fermionic modes. 

It is helpful to check the corresponding scaling behavior as a function of temperature for the first 
few thermodynamic potentials. The internal energy of the canonical ensemble, U = — (dW/d(3) v , 
displays the same scaling behavior as the free energy. The entropy is given by the expression: 



S = -13 2 




= -(3 2 



(3- 2 W((3) - (3 



'dW s 



(53) 



Finally, since the specific heat at constant 



(54) 



and we infer that it scales at high temperatures as (3 1 
volume is given by: 

*--'(£), 

we infer that it also scales as (3 1 at high temperatures. The corresponding scaling behavior at 
low temperatures, in agreement with the expected result for the low energy gauge theory limit 
can be extracted, as in the previous section, by considering instead the t^oo asymptotics of these 
expressions. 
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5 Conclusions 



We have shown that all of the perturbatively renormalizable and anomaly-free supersymmetric string 
theories with Yang-Mills gauge fields: heterotic, type I, and type I' admit stable and tachyon-free 
finite temperature ground states in which we can formulate an equilibrium statistical mechanics of 
strings in the weak coupling limit. Preceeding attempts [33, 26, 40, 36, 41, 43, 42] to formulate an 
equilibrium statistical mechanics with a tachyon-free supersymmetric string canonical ensemble have 
failed for a variety of reasons; a detailed account appears in section 2 of hep-th/0409301vl. These 
works do not correctly incorporate the Euclidean T-duality transformations linking the thermal 
vacua of the supersymmetric string theories in pairs. Neither do they meet the infrared consistency 
conditions we have required in our analysis, matching self-consistently with the low energy finite 
temperature field theoretic physics. In addition, string theoretic consistency conditions, such as 
modular invariance, have often been violated. The results in this paper demonstrate irrefutably 
that the widespread misconception that the canonical ensemble of superstrings necessarily breaks 
down beyond a limiting (Hagedorn) temperature is simply wrong. Finally, in the case of the 
pure type IIA and type IIB thermal ground states, we show that it is not possible to eliminate 
low temperature tachyons from the thermal spectrum while meeting the constraints from modular 
invariance. This fact also precludes interpretation of thermal winding mode tachyons as indication 
of a string scale "Hagedorn" phase transition in the type II theories: the canonical ensemble cannot 
be defined in the full temperature range, starting at T=0. 

For each of the six superstring ensembles, we have shown that the growth of the free energy at 
high temperatures far above the string scale is only as fast as that in a 2d quantum field theory: 
F(T) ~ T 2 . This behavior exemplifies the equivalence of perturbatively renormalizable superstring 
theories to the 2d gauge theory of diffeomorphism and Weyl invariances. We have simultaneously 
verified that the low energy field theoretic limit of our expressions for the string free energy in 
each case recovers the expected T 10 growth characteristic of a 10D quantum field theory, at low 
temperatures where only the lowest-lying massless field theory modes contribute. These results 
follow directly from an explicit evaluation of the one-loop modular integrals. The T 2 growth of the 
high temperature string free energy was originally conjectured by Atick and Witten in [41]. We 
have shown that the T 2 scaling can also be inferred from the thermal duality relations linking the 
thermal ground states of the superstring ensembles in pairs. This is in precise analogy with the T 2 
scaling of the free energy of the closed bosonic string as a consequence of its thermal self-duality, 
as shown by Polchinski in [1]. In fact, we can infer the existence of a duality phase transition in the 
Kosterlitz-Thouless universality class mapping the finite temperature ground state of the E 8 xE 8 
heterotic string to its Spin(32)/Z 2 Euclidean T-dual, analogous to the self-duality phase transition 
observed in the closed bosonic string ensemble [20]. The Kosterlitz-Thouless universality class is 
characterized by an infinite hierarchy of thermodynamic potentials displaying analyticity at the 
critical temperature, in this case, T c = 1/2W 1 / 2 [18, 35, 34, 20]. 

We will close with mention of an important insight that applies more broadly to the development 
of a fundamental, and nonperturbative, formulation for String/M Theory. We remind the reader 
that perturbative string theory as formulated in the worldsheet formalism is inherently background 
dependent: the "heat-bath" representing the embedding target space of fixed spatial volume and 
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fixed inverse temperature is forced upon us, together with any external background fields char- 
acterizing the target spacetime geometry. Thus, we are ordinarily restricted to the canonical en- 
semble of statistical mechanics. We should caution the reader that while an immense, and largely 
conjectural, literature exists on proposals for microcanonical ensembles of weakly-coupled strings 
[25, 40, 31, 42, 43], the conceptual basis of these treatments is full of holes. Some of the pitfalls 
have been described in [14, 41]. One could argue that, strictly speaking, the microcanonical en- 
semble is what is called for when discussing quantum cosmology, or the statistical mechanics of 
the Universe [49]: the Universe is, by definition, an isolated closed system, and it is meaningless 
to invoke the canonical ensemble of the "fundamental" degrees of freedom. However, there remain 
many simpler questions in both early Universe cosmology, and in finite temperature gauge theories, 
that are approachable within the framework of the canonical ensemble. Limited use of the micro- 
canonical ensemble under certain assumptions is also possible, but the constraints imposed by the 
Jeans instability, and by thermal back reaction, need to be kept in mind. 
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